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Abstract 

We construct static and time dependent exact soliton solutions for a theory 
of scalar fields taking values on a wide class of two dimensional target spaces, 
and defined on the four dimensional space-time S"^ x IR. The construction is 
based on an ansatz built out of special coordinates on S"^. The requirement for 
finite energy introduces boundary conditions that determine an infinite discrete 
spectrum of frequencies for the oscillating solutions. For the case where the target 
space is the sphere S'^, we obtain static soliton solutions with non-trivial Hopf 
topological charges. In addition, such hopfions can oscillate in time, preserving 
their topological Hopf charge, with any of the frequencies belonging to that 
infinite discrete spectrum. 



1 Introduction 

In this paper we consider a non- linear theory of a complex scalar field u, on the space- 
time 5*^ X IR. Our considerations apply to a wide class of target spaces, but the case 
of most interest is that of the sphere S"^, in which case the field u parametrizes a plane 
corresponding to the stereographic projection of S"^. The action is the integral of the 
square of the pull-back of the area form on the target space. Therefore, it is quartic in 
derivatives, but only quadratic in time derivatives. The theory is integrable in the sense 
that it possesses a generalized zero curvature representation of its equations of motion 
and an infinite number of local conservation laws ^1 12] . The conserved currents are 
associated to the invariance of the theory under the area preserving diffeomorphisms 
on the target space. 

We construct an infinite number of static and time dependent exact soliton solutions 
using an ansatz j21 12] that reduces the four dimensional non-linear equations of motion 
into linear ordinary differential equations for the profile function. In the case of the 
target space S'^ the solitons carry non-trivial Hopf topological charges. Although the 
topology is the same as that of other models possessing hopfion solutions |S1 IH El, the 
Derrick's scaling arguments is circumvented in a different manner. The stability of the 
static three dimensional solutions comes from the fact that the physical space is S^ and 
that introduces a length scale given by its radius. A model in Euclidean space where 
a similar stability mechanism occurs is discussed in jH]. 

The requirement for finite energy leads to boundary conditions that determine an 
infinite discrete set of allowed frequencies for the oscillations of the solutions. Those 
solutions can be linearly superposed since the profile function satisfies a linear equation. 
It turns out that the energy of the superposed solution is the sum of the energies of 
the modes, and in this sense the modes are decoupled. However, the profile function 
can take values only on some intervals on the real line, which depend on the choice 
of the target space. Therefore, not all superpositions are allowed and that introduces 
some sort of coupling among the modes. The allowed superpositions for the case of 
the target space S^ are discussed in detail. One of the most interesting superpositions 
corresponds to oscillating hopfion solutions. We show that it is possible to superpose to 
the static hopfion soliton any number of oscillating modes with any of the frequencies 
belonging to that infinite discrete spectrum. Although the solution oscillates in time, 
its topological Hopf charge is preserved. The only constraint on such superposition 
appears on the intensity of each mode. 

The paper is organized as follows: in section |21 we introduce the model and discuss 
its integrability properties, in section 01 we propose the ansatz and construct the exact 
soliton solutions, in section 0] we discuss the energy, boundary conditions and allowed 
frequencies. The case of the target space S'^ is discussed in section El 



2 The model 

The metric on the space-time S*^ x ]R is given by 

ds^ = df-rl i———- + {l-z)dipi^ + zd^2^j (1) 

where t is the time, z, and cpi, i = 1,2, are coordinates on the sphere S^, and < 2; < 1, 
< V'i < 27r, and tq is the radius of the sphere S^. Embedding S^ on IR"^ we get that 
the Cartesian coordinates of the points of S^ are 



Xl = roy/zCOSip2 X3 = Tq y/l — Z COS (yJi 



X2 = rQ^/zsm^p2 x^ = Tq \/1 — z sinipi (2) 

The model is defined by the action 



h 



S= dt dT.^ (3) 

where the volume element on S^ is dH = ^ dzdipi dip2, and 

hf,u = d^ud^u* - dyud^u* (4) 

where d^ denotes partial derivatives with respect to the four coordinates on S^ x IR, 
namely C^ = (^; ^? V'l; ¥^2), m is a complex scalar field, and 7 = 7(| -u p), is a real 
functional of the squared modulus of u, and it defines the geometry of the target space. 
The metric on target space is given by 

, r, dudu* , , 

da^ = 5 

7 

Some cases of interest are the following: a j 7 = 1 corresponding to the plane with 
coordinates being the real and imaginary parts of u, b) 'y = (1— | m P) , with | u p< 1, 
corresponding to the Poincare hyperbolic disc, and cj 7 = (1+ | u p) corresponding 
to the sphere 5*^. In such case, u is related to the three dimensional unit vector n 
{n^ = 1) defining the sphere, through the stereographic projection 

n = ■ (u + u*, -i iu- u*) ,\u P -1 ) (6) 

1+ I M p ^ ^ 

Notice that in the case of the target space being S"^ , the action Q corresponds to the 
quartic term of the Skyrme-Faddeev action jl] which has been studied on the space-time 
53 X IRin [H. 

The Euler-Lagrange equations following from Q are given by 



together with its complex conjugate, and where 

K^ = h^i^d^u = (dudu*) d^u — {du) d^u* (8) 

The model Q possesses an infinite number of local conserved currents given by 

K,6G K; 6G 

7 OM 7 ou* 

where G is a function of u and u* but not of its derivatives. Using the equations of 
motion ((7j) and the identities 

K^ d^u = K^ dV - Kl d^u = (10) 

one can check that is indeed conserved, i.e. d^^Jf^ = 0. The symmetries associated to 
such conservation laws are the area preserving diffeomorphisms of the target manifold. 
Indeed, the tensor h^y/'~^ is the pull-back of the area form on the target space 

dA = t (11) 

7 

Therefore, the action (jHl) is invariant under diffeomorphisms preserving the area (jllll . 
and (P) are the corresponding Noether currents jHl 12] ■ 



3 The ansatz 

Following [21 El we introduce the ansatz 

u = f{t, z) e^('"i^i+'"2^2) (12) 

where rrii, i = 1, 2, are arbitrary integers, and / is a real profile function. Replacing it 
into the equation of motion Q one gets 



1 ) rlVt \ 1 



dt{-^\--^^^d4n^-^\=Q (13) 



where 

VL = mlz + mlil - z) (14) 

We now make a change of variable in the profile function, introducing a function g by 

dp 
dgr-.^ (15) 

7 

For instance, in the particular cases discussed below (0) we get that: a) for the plane 
where 7 = 1 one has g = f^, and g > 0; b) for the Poincare hyperbolic disk where 
7 = (1— I M p) , one has g = 1/ (1 — /^), and g > 1', c) for the sphere 5^ where 
7 = (1+ I u p) one has 5- = 1/ (1 + p), and < 5- < 1. 



So, with the change (fT3j) one gets that (fT^ becomes a hnear partial differential 
equation in g 



ro^fi 







The solutions can be obtained by separation of variables introducing 

g{t,z) = J{t)H{z) 
and then J and H have to satisfy 



and 



z{i-z)d.^ 





dlJ + u 


'J = 


[{q'z + {l- 


-z))d,H 


^2, ,2 

4 



[q^z + {l-z)) H = 



where uj"^ is the separation of variables constant, and where q is 



Q = 



mi 
m2 



(16) 
(17) 

(18) 
(19) 

(20) 



The equation (fTUj) is what is called a Heun equation 0. It is a generalization of 
Gauss hypergeometric equation in the sense that it has one extra regular singular point. 
Inddeed, flT^ reduces to the hypergeometric equation for g = 1. Notice that (fT^ is 
invariant under the joint transformations q <-*> 1/g and z ^^ 1 — z. Therefore, if H„„j (z) 
is a solution of (fTT?|l for some value of q, so is 



^q^LO 



Hl/q,ui {z) = Hq^^ (1 - Z) 



(21) 



for the inverse value 1/g. One can obtain solutions of ()19|1 in powers series around 
z = Q and for < g < 1. The solutions for g > 1, are then obtained from those using 
the above symmetry. The power series solutions are given by 



Hq,Lo {z) 



-^q,LU 



I oo 

\ ra=2 



(22) 



where the positive constants Vq^^ are chosen such that the maximum absolute value of 
^q,uj {.z) in the interval < z < 1, is unity. The coefficients, for n > 2, are determined 
by the recursion relation (with cq = 0, and Ci = 1) 



n (n 



2 ^2 



00 r, 







4 



[n 



-71-2 



(23) 



+ 



'uj'^r? 



0_(^_l)(^_2) + (g2_l)(^_l)2 ^^^_ 



We will be interested in solutions with u real, and therefore the solutions of ()18j) are 
trigonometric functions. Of course, if ^f is a solution of (HD), so is a g + P, with a and 
/3 constants. We will then use the following normalization for the solutions of ()16|) 



9q,i 



[sin (cut + 5) Hg^^{z) + 1] 



cu^O 



(24) 



with Hq^^ {z) given by ()22|) . The advantage of such normahzation is that the solutions 
(|24p take values on the real line from zero to unity only. Therefore, they are admissible 
solutions for the case where the target space is the sphere S'^ (see discussion below 

We normalize the static solutions of ()16|) as 

_ \Ti{q^z + l-z) 
9,,o j^, (25) 

So, Qqfi is a monotonic function varying from zero at z = to unity aX z = 1. Notice 
that Qqfi -^ z as q ^ 1. A decreasing function can be obtained by the interchange 
gq,o — ^ 1 — gqfl- Therefore, (j^ are also admissible solutions for the case where the 
target space is the sphere S*^. 

The admissible solutions for the cases of the plane or the Poincare hyperbolic disk 
can then be written as 

„plane ^ ^ ^Poinc.disk ^ ^ i -i fna\ 

9q,^ = " 9q,u 9g^^ = a Qg^^ + 1 (26) 

with a being a real and positive constant, and Qg^^^ being given either by (j^ or (J2S)). 

4 The energy and boundary conditions 

The energy for the solutions obtained through the ansatz (fT^ . p5|) . is given by 

E= f d^n = 16n'ro f dz Q ( ^^f + \ (d^g)'] (27) 

JS^ Jo \4:Z (1 — Z) Tq J 

where Ti is the Hamiltonian density associated to Q- Using the equation of motion 
(fTB|) one gets that the energy of static configurations is 

^static = ^ {ngd,g)\lzl (28) 

Therefore, the energy of the static solutions ()25|1 is 

„ / N 16 7r^ , , (q — l/q) , ^ 

^(^.,o) = Imimal ^ , '^^ (29) 

Notice that the solutions of (fT^ . for cu 7^ 0, which do not vanish at z = or z = 1, 
have a logarithmic divergence on its first derivative, at those points. We then observe 
from (jTfjl that such solutions do not have finite energy. Consequently, we shall impose 
the following boundary conditions 

i/g,, (0) = /7,,^ (1) = for ^7^0 (30) 



In addition, using the equation of motion ()16|) one obtains that 

dE 327r2 



dt tq 



(ndtgd.g) llzl, (31) 



Therefore, the energy is conserved for the static solutions, and for those solutions (j24j) 
satisfying the boundary conditions (jSOI)- 

Multiplying (fT!I|l by Hq,j and integrating in z, one gets that 

o'^"4707)^'-"^r^"^ Kj' = -4 {^H,,^K^) \t=o (32) 

Consequently, the energy ^I7\ for the solutions (j^ satisfying ^^ is 

^, , 16 7r^ . . . , , 

E [gq,uj) = I mi m2 I A,,,^; (33) 

where 

(34) 
The fact that A^^^^ is invariant under q ^^ 1/q, follows from the symmetry (PT|) . In Table 
121 we give the values of A^^^ for some chosen values of q and some allowed frequencies 
u. So, the energies (^^ and (jHSl) do not depend upon the signs of the integers mi and 
1712 and it is invariant under the interchange mi ^^ m,2 (see (|2Up). 

Multiplying ()19|1 by Hq^Q, subtracting from the same relation with u interchanged 
with uj, and integrating in z, one gets that 



2 -2 

u — u 



Jo ^^ Az (1 - z) ^^'"^ ^'^•'' " ^ (^ (^'^•'^ ^^'"^ " ^^'^ ^^''')) '^=° ^^^"^ 



Similarly, differentiating ()19|1 w.r.t. z once, multiplying by -f^^j;^, subtracting from the 
same relation with u interchanged with lD, and integrating in z, one gets that 

{uj'-iu') fdznH'q^^H'q^^ = \{Az{l-z) n {H'q^^ H'l^ - H'l^ H'q^^)) l^zj (36) 

JO Tg 

Consequently, for the solutions P^ satisfying (jHUI) . and having finite first and second 
2;-derivatives at z = and z = 1, one gets the orthogonality relations 

I' dz ^^ ^_ ^^ Hq,^ Hq,^ = f^ dz Q g;^ i/;^ = for u^u (37) 

It then follows that the energy (j?r|) of a linear combination of solutions P^ satisfying 
finn|) . is just the sum of the energies of each solution, i.e. 

E {A Qq,^ + B gq,^) = A^E (gq^^) + B^E (gq^,,) for lj ^ cD (38) 



Therefore, in that sense, the modes for the same value of q are decoupled. However, the 
intensities in which they enter in the superposition are not independent. As discussed 
below (J15p the real values that the profile function g can take depend on the target 
space under consideration. So, when we take linear combinations of the solutions we 
have to respect those constraints. In section El we discuss those constrained linear 
combinations in detail in the case where the target space is the sphere S"^. 

The boundary conditions ()30|) lead to a discrete spectrum of allowed frequencies uj. 
Indeed, for the case where g = 1 the series (j^^ truncates whenever 

^^ = n(n+l) n = 1,2,3,... (39) 

and the corresponding polynomials satisfy the boundary conditions (jHOJ. The first four 
of those polynomials are given by (with the normalization as in ()22|) ) 

H^Q = gVs [z - 3z^ + 2z^) (40) 

Hi^i2 = 16 [z - Gz"^ + lOz^ - 5z^) 

2450 {z - lOz^ + 30z^ - 35z^ + Uz^) 



Hi 



20 



VT (3^5 + 5^6) V15 - 2^30 

where the first index refers to g = 1 and the second to n{n + 1) as given by ()39|). 

For g 7^ 1 the series (j^ does not truncate, and the frequencies u leading to 
solutions satisfying ()30p can be found numerically. We give in Table ^ the first three 
frequencies for some chosen values of q. In addition, the same frequencies hold true 
under the exchange q ^^ 1/q due to the symmetry (J2H). Therefore, the frequencies 
become smaller as q departs from unity, either to smaller or greater values than unity. 

In Figure n we exemplify the shape of the functions Hq^^, by plotting the first three 
modes Hg^^^, for q = 3/4 and q = 1/10 and the frequencies Ui, i = 1,2, 3, given in Table 
UJ Due to the symmetry (PT|) . Hi/g^^. can be obtained by reflecting the plots around 
z = 1/2. Notice that the polynomial solutions for q = 1, given in (jlUI) . are invariant 
under fETj) . One observes that as q decreases, the functions Hg^^ get deformed in a way 
that their first derivatives at z = 1 increase. For g > 1, it follows, due to (J2H) . that the 
first derivative of Hg^^ at ^ = increases as q increases. 

5 The case of S"^ as target space 

According to the comments below ()15|) . in the case where the target space is the two 
dimensional sphere S^, we have that 7 = (1+ | m P) , and so g = 1/ (1 + p) and 
< g < 1. The ansatz ()12|) becomes 



u = J- — - e^('"i^i+™2V2) (41) 

7 



As we have discussed, the solutions (j^ and (^3)) are admissible solutions for the 
case where the target space is S"^. However, one can construct more admissible solu- 
tions, i.e. with < g < 1, hj taking linear combinations of the solutions (|24j) and 
fl25|) . An interesting case corresponds to linear combinations of the static solution Qq^o 
given by ()25|) . with one or more time dependent solutions of the type (j23)- It leads to 
oscillating (in time) hopfion solutions. Let us consider the case g < 1 first. The static 
solution ggfi vanishes at z = 0, and has z-derivative equals to (g^ — 1) / Ing^ there. The 
function fg (z) = [(g^ — 1) /Ing^] Vg^^Hg^^ (z), with Hg^^ given by (j^ and satisfying 
(IHUj) . has the same behaviour at ^ = as gg^. In addition, C2 given by (J2S1), is negative 
for the case of Hg^^ satifying (jHUj) . and so fg{z) grows slower than ggQ for small z. 
We do not have a rigorous proof, but by careful direct inspection we found that, as z 
increases from zero to unity the absolute value of fg (z) never exceeds the value of gg^, 
which is a positive monotonic function of z. It then follows that the combination 

(g2 _ 1) 
gg,o + a— i — 5— Vg^u; sm {ut + 6) Hg^^ (z) (42) 

with < a < 1, takes values between zero and unity only, and so it is an admissible 
time dependent solution for the target space S'^. Therefore, by adding up such of type 
of solutions and dividing by their number, one gets that the time dependent solutions 

f,^^ , , In (q"^ z + 1 — z) 1 (g^ — 1) .^-^ . , r ^ ^^ / x / ,„x 

4^ (^' ^) = i^^2 — - + iv 1^ ^ "^ '''''" '''' ^^'^ + ^'^ ^^'-' ^'^ ^^^^ 

where N is the number of modes (frequencies) entering into the sum, and a^'s are real 
coefficients satisfying < Oj < 1, are admissible solutions for the target space S'^. The 
energy of such solutions, according to (j^ . (j^ . (j^ . (jH^ and (jHHj) . is given by 

with Ag,<^^ given by dSH). 

Solutions of the type (jlHj) for g > 1 can be obtained using the symmetry ((211), 
i.e. g[i^ it,z) = g^^^ (t, 1 — z). Solutions that decrease from unity at 2; = to zero 
at 2; = 1 can be obtained by the symmetry g^^^ it,z) ^ 1 — g^^^ it,z). The energy 
fl44|) is invariant under those two type of transformations. As we now show, all these 
solutions correspond to oscillating hopfions, i.e. solutions that oscillates in time and 
have a constant (in time) non-trivial Hopf number. 

For any fixed time t our solutions define a map from the physical space S^ to the 
target space 5*^, and so it is a Hopf map ^^. We now show that the Hopf invariant (the 
linking number) is independent of time for the admissible solutions we have constructed. 
In order to calculate de Hopf index of the solution we introduce 



*i = ^/g cos{'m2ip2) ^3 = Jl- g cos{miipi) 



*2 = -^sin(m2V?2) ^4 = Jl - g sm{miipi) (45) 



which defines another 3-sphere Si, since $^ + $| + $^ + $^ = 1. The field u in (jiTj) 
can be written as 

^ = -^ ^^r (46) 

Since u parametrizes the sphere S"^ through the stereographic projection (jH)), we have 
that (jinjl gives the map S^ -^ S"^. So, the Hopf index is in fact evaluated through the 
map S^ ^ S^ ^ S"^, as we now explain. We introduce the potential 

A = -(^Z^VZ-VZ^ Z^ (47) 



where 



Z=l^']- Zi = $3 + «$4 ; ^2 = $1 + i$2 (4J 



^2 



and the differential operator V is the gradient on the physical space S^. The Hopf 
index is defined by the integral ^U] 

Q^= I dT. A- curl A (49) 

where dH is the volume element on the physical space S^. Evaluating we get 

7 {1 — g) ^ g ^ 

A = -mi — = e^^ +1712^= e<^2 (50) 

and 

curM = 2 (9^5f (-m2A/r^^e^i +miV^e^2J (51) 

Consequently we get that 

QH = mm2[9{t,z = l)-g{t,z = 0)] (52) 

Therefore, the solutions (j^^ have Hopf index 

Qh {gqfi) = mim2 (53) 

and so they are static hopfion soliton solutions. The time dependent solutions (j^^ 
are constant at z = and z = 1 due to the boundary condition (jSH)), i-e. gqu (^,0) = 
9qu] (^) 1) = 1/2. Therefore, they carry no Hopf number 

Qh {gq,u.) = o; 7^ (54) 

The solutions (j^Hj) although time dependent, also have contant values at 2; = and 
z = 1, determined by their static component. It then follows that 

Qh {gi""^) = 77117712 (55) 

and so they do correspond to oscillating hopfion soliton solutions. 



The Hopf index can also be calculated as the linking number of the pre-images 
of two points of S'^ [TU]. Notice, from (jH)), that the north pole of S"^, n = (0,0,1), 
corresponds to m ^ cx), and so from (j4ip to g = 0. On the other hand, the south pole 
of S*^, fi = (0, 0,-1), corresponds to u = 0, and so to (? = 1. Therefore, from (J13j) and 
PH|) . we see that the pre- image on S^ of the north pole of 5^ correponds to Zi = e^"^^'^^ 
and Z2 = 0, whilst the pre-image of the south pole to Zi = and Z2 = e"*™'^'^^. For 
the static solution (j^ and the time dependent solutions (ji^ . the pre-images in the 
spacial 5*^ of these two circles in 5*1 are constant in time. In addition, those two circles 
in S"! pass through each other mim2 times as (pi and (^2 varies from to 2n in S^. 
So their linking number is mim2, and that is the Hopf index. For the time dependent 
solutions ()24|) it is not possible to have g = and g = 1 on different points on the 
spatial S*^ at the same time t. Therefore, the pre-images of the north and south poles 
of S"^ never link, and so the Hopf index of such solutions vanishes. 
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Q 


u;lrl/A 


ujlrl/A 


ulrl/A 


1 


2 


6 


12 


3/4 


1.983765266031 


5.987866563654 


11.988820961965 


1/2 


1.913643550920 


5.923271548463 


11.930495946441 


1/4 


1.733644829967 


5.672599570202 


11.668111143086 


1/10 


1.517649738276 


5.238029617258 


11.066608078277 



Table 1: The first three frequencies leading to solutions (j^^ satisfying the boundary 
conditions (jHUj) for some chosen values of q. Due to the symmetry ()21|) the frequencies 
are invariant under the interchange q ^^ 1/g. 



(1 


^^q,uii 


^^q,u)2 


^^q,u-i 


1 


4/3 


27/5 


64/7 


3/4 


1.362643 


4.693750 


9.501984 


1/2 


1.505658 


4.275566 


9.428440 


1/4 


2.065180 


4.819658 


8.855009 


1/10 


3.739188 


7.930185 


13.004162 



Table 2: Numerical values of A^^^-, as defined in (jS^), for some chosen values of q and 
for the three frequencies given in Table [T] Due to (PT|) , A^ ,^ is invariant under the 
interchange q ^^ 1/q. 



q 


■^g,!^! 


Vq,iLj2 


Vq,ui3 


1 


1/4 


i/eVs 


1/16 


3/4 


0.28518993 


0.11915683 


0.07078933 


1/2 


0.33261165 


0.15277952 


0.08701028 


1/4 


0.40999103 


0.20414992 


0.12695732 


1/10 


0.50274925 


0.25672647 


0.16731092 



Table 3: Numerical values of the normalization constant Vq^^Ji, as defined in ((221), for 
some chosen values of q and for the three frequencies given in Table ^ Due to ()21|1 , 



Vq^uj is invariant under the interchange q 



1/q. 
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q=3/4 



0.5 



-0.5 



0.5- 



-0.5- 




Figure 1: Plots of if, 



q^uji\z), i = 1,2,3, as normalized in (j221), for 
q = 3/4 (top) and q = 1/10 (bottom) and for the three 
frequencies given in Table H The number of zeroes increase 
with the increase of cuj. The plots of Hg,^ and Hi/g^^ are 
related by reflection around z = 1/2 due to the symmetry 



12 



References 

[1] O. Alvarez, L. A. Ferreira and J. Sanchez Guillen, "A new approach to integrable 
theories in any dimension," Nucl. Phys. B 529, 689 (1998) ||arXiv:hep-th/971014f| . 

[2] O. Babelon and L. A. Ferreira, "Integrability and conformal symmetry in higher 
dimensions: A model with exact Hopfion solutions," JHEP 0211, 020 (2002) 
|arXiv:hep-th/0210154j . 



[3] H. Aratyn, L. A. Ferreira and A. H. Zimerman, "Exact static soliton solutions of 
3+1 dimensional integrable theory with nonzero Hopf numbers," Phys. Rev. Lett. 
83, 1723 (1999) |arXiv:hep-th/9905079|; and "Toroidal solitons in 3+1 dimen- 
sional integrable theories," Phys. Lett. B 456, 162 (1999), arXiv:hep-th/990214'l] 

[4] L. Faddeev, in 40 Years in Mathematical Physics , (World Scientific, 1995); 
L. D. Faddeev and A. J. Niemi, "Knots and particles," Nature 387, 58 (1997) 
|arXiv:hep-th/9610193j ; R. A. Battye and P. M. Sutcliffe, "Knots as stable soliton 
solutions in a three-dimensional classical field theory," Phys. Rev. Lett. 81, 4798 
(1998) |arX iv:hep-t h/9808129 ; J. Hietarinta and P. Salo, "Ground State In The 
Faddeev-Skyrme Model," Phys. Rev. D 62, 081701 (2000). 

[5] D. A. Nicole, "Solitons With Nonvanishing Hopf Index," J. Phys. G 4, 1363 (1978). 

[6] L. A. Ferreira, "Euclidean 4d exact solitons in a Skyrme type model", 
|arXiv:hi^th/0406227 | 

[7] R. S. Ward, "Hopf solitons on S(3) and R(3)," Nonlinearity 12 (2) (1999) 241-246, 
|a7Xiv:hep-th/9811 176.1 . 

[8] L. A. Ferreira and A. V. Razumov, "Hopf solitons and area preserving diffeomor- 
phisms of the sphere," Lett. Math. Phys. 55, 143 (2001) ^arXiv:hep-th/0012176| . 

[9] A. Ronveaux, Heun's Differential Equations, Oxford University Press (1995) 
S. Yu. Slavyanov and W. Lay, The Heun Class of Equations, chap. 3 in Special 
Functions: A Unified Theory Based on Singularities, pp. 97-162, Oxford University 
Press (2000). 

[10] R. Bott and L.W. Tu, Differential Forms in Algebraic Topology, Graduate Texts 
in Mathematics 82, Springer- Verlag (1982). 



13 



